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At the University of Arizona, we have a major effort to study core-collapse supernovae, 
neutrino opacities and transfer, multi-dimensional radiation hydrodynamics, pulsar kicks, and 
gravitational radiation from supernovae. We have developed a precision ID (spherical) multi- 
k*v group implicit Boltzmann solver, SESAME 1 [1, 2], using the Feautrier technique, with the 

tangent-ray method to achieve excellent angular resolution. This code provides the spherical 
benchmark and testbed for our multi-D efforts (e.g., VULCAN/2D, see below). 

In an effort to explore non-standard effects quickly and heuristically, Thompson, Quataert, 
and Burrows [3] (TQB) have modified SESAME to incorporate the centrifugal effects of dif- 
\ ferential rotation in an approximate way, as well as the heating and angular momentum 

redistribution effects of viscosity (perhaps due to the magneto-rotational instability, MRI). 
£^ ■ The neutrino-driven mechanism is already close to explosion and requires an increase in the 

energy deposition rate in the "gain region" of only about 25% to 40%. TQB have shown that 
rapid rotation of the initial Chandrasekhar core can provide such a boost through the agen- 
cies of centrifugal support and viscous heating alone. As the explosion commences, residual 
neutrino cooling subsides and neutrino heating takes over, but the initial boost due to viscous 
heating in these calculations helped trigger explosion. They have also shown that a delayed 
Qh' explosion more naturally ejects high Y e material (~ 0.5) than a prompt (never achieved in 

realistic simulations) or early delayed explosion. Furthermore, the larger residual core mass 
from which the neutron star is formed is more in keeping with larger measured pulsar masses. 
However, a core that explodes using a viscous boost is also rotating more rapidly than the in- 
ferred birth spin rates of known pulsars. Hence, a mechanism that sheds angular momentum 
during and after explosion is required. TQB show that here too the viscosity that heats the 
core can also lead to the diffusion of angular momentum out of the core (by a factor of ~5 
by the end of their calculations) on a timescale of seconds. These calculations, being spheri- 
cal, are by no means definitive, and true multi-D simulations, including magnetic fields, are 
required to answer the questions the work of TQB pose. However, these results are intriguing. 

As a biproduct, TQB have also estimated the evolution of a core's rotational profile 
during collapse and after bounce, i.e. the mapping between the initial and the "final" angular 
velocities. This study with neutrino transport complements that of Ott et al.[4] without 
neutrino transport (which nevertheless explored a wide range of initial conditions). 

As TQB show, a useful criterion for explosion is that the timescale (r a ^) for matter to 
flow from the stalled shock into the cooling region near the neutrinosphere is larger than 
the characteristic heating timescale (tjj)- This criterion is loosely equivalent to the L u vs. 
M criterion of Burrows & Goshy[5]. For a given M, there is a threshold L u , above which 
explosion ensues. That threshold is lower in 2D/3D than in ID. In ID, r n ^ is too small and 
matter settles onto the core before sufficient heat is deposited by the neutrinos. However, as 
shown by Burrows, Hayes, &; Fryxell[6] (BHF), Herant et al.[7], and Janka and Miiller[8], in 
2D with the neutrino-driven overturning instability and convective motions behind the stalled 
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shock enabled, T ac [ v can be larger than th- The swirling motions keep the matter in the gain 
region longer, increasing T a ^ v . As a result, the 2D calculations of BHF and Herant et al., 
without rotation, led to explosions. Note that if matter did not settle onto the core at all, or 
if the cooling rates near the neutrinospheres interior to the gain region were low, explosion 
by neutrino heating would naturally, and trivially, result. We know how to obtain explosions. 
The question is whether Nature agrees. 

Importantly, the results of BHF and Herant et al. were obtained using simple neutrino 
transfer and transport algorithms. Though they may have been substantially correct, it is 
thought that better calculations of the neutrino-matter coupling, particularly in the semi- 
transparent gain region, are required to falsify and/or check the convective neutrino-driven 
mechanism of explosion that I, for one, think obtains. Hence, a number of groups (for instance, 
at Arizona, MPA, and ORNL) are developing more sophisticated multi-dimensional, multi- 
group neutrino transfer schemes and in the next few years we should be able to make a much 
more detailed assessment of this mechanism. 

In this spirit, we [9] have developed the code "VULCAN/2D" and are testing it in full 
Boltzmann and MGFLD modes. Figure 1 provides a snapshot in a 2D MGFLD run of the 
electron neutrino flux vector field at 7.8 MeV, superposed on an entropy color map. VUL- 
CAN/20 is an ALE code with remapping that implicitly solves the 6-dimensional problem 
(l(time) + 2(space) + 2(angles) + 1 (energy-group)) using the S n method to tile angular space. 
It is parallelized using MPI, is flux-conserving, and smoothly matches to the diffusion regime, 
but does not include the Doppler and aberration velocity terms. The former is important, 
particularly around bounce and will be included explicitly in operator-split fashion in a later 
version. However, the achievement of the current VULCAN/2D is a milestone, being the first 
such time-dependent code available in core-collapse studies with truly "2D", as opposed to 
ray- by-ray or only spherical, capability. Ray-by-ray codes, such as used by BHF, Rampp & 
JankaflO], and Buras et al.[ll], can exaggerate the variation in the neutrino energy density 
and flux with angle by as much as 20%-50%, with the concommitent artificial variation in the 
heating rates with angle. However, it is not yet known how, in the integral sense of global 
heat deposition, this difference might effect the mechanism itself. 

Whether rotation is central to the mechanism of core-collapse supernovae is not known. As 
the work of BHF, Herant et al.[7], and Janka and Miiller[8] suggest, it may not be. However, 
as TQB, Ott et al.[4], and Fryer & Heger[12] demonstrate, rotation (naturally amplified during 
collapse and compaction) introduces interesting new phenomena. In addition to decreasing 
the effective gravity against which the explosion must emerge, enlarging the gain region, and 
generating funnels along the poles through which an explosion can more easily pierce [13], 
stars do in fact rotate. Unfortunately, how much the Chandrasekhar core rotates just before 
collapse is currently a mystery. An example of such post-bounce funnel structures is given 
in Fig. 2. (I refer the reader to the bi-polar structure that HST has revealed in SN1987A.) 
Therefore, the consequences of the fact of rotation, and the possible synergistic effects with 
magnetic fields[3, 14], deserve further scrutiny in the context of the neutrino-driven paradigm 
of core-collapse supernova explosions. 
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Figure 1: A snapshot 5.6 milliseconds after bounce. The 2D MGFLD variant of VULCAN/2D 
was used. The color map is of the entropy and the vectors are the electron-neutrino differential 
fluxes at 7.8 MeV. The core is not rotating, but in the beginning of the calculation it was 
given a modest anisotropic perturbation in the density at the north pole. The approximate 
scale is 300 kilometers. (Obtained in collaboration with R. Walder.) 
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